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(57) Abstract: A method of determining parameter values used in 
impedance analysis of a subject- The method includes using a pro- 
cessing system to determine a number of impedance measurements 
at a corresponding number of frequencies. The impedance measure- 
ments are used to determine estimates of parameter values, with the 
estimates being used to determine theoretical impedance values based 
on the parameter value estimates. The theoretical impedance values 
are compared to the measured impedance values to allow a modifica- 
tion direction to be determined for at least one of the parameter value 
estimates. This is then used to modify at least one parameter value 
estimate, with the parameter values being determined at least in part 
from one or more modified parameter value estimates. 
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IMPEDANCE PARAMETER VALUES 
Background of the Invention 

The present invention relates to a method and apparatus for determining parameter values used in 
impedance analysis of a subject 

5 Description of the Prior Art 

The reference in this specification to any prior publication (or information derived from it), or to any 
matter which is known, is not, and should not be taken as an acknowledgment or admission or any 
form of suggestion that the prior publication (or information derived from it) or known matter forms 
part of the common general knowledge in the field of endeavour to which this specification relates. 

10 One existing technique for determining biological parameters relating to a subject, such as cardiac 
function, or oedema, involves the use of bioelectrical impedance. This involves measuring the 
electrical impedance of a subject's body using a series of electrodes placed on the skin surface. 
Changes or values of electrical impedance can then used to determine parameters, such as changes in 
fluid levels, which can in turn be used as indicators of attributes of cardiac cycle or oedema. 

15 Complex signal processing is required to ensure measurements can be interpreted. An example of 
this is described in Intemational patent publication no WO2004/032738, in which the responsiveness 
of a patient to an applied current is measured to determine a number of impedance values. The 
measured impedance values are plotted as a frequency dependent locus, which is then extrapolated to 
determine impedance values at zero and infinite applied frequencies. 

20 However, such plotting and extrapolation techniques are computationally expensive, making it 
difficult to analyse impedance measurements in real time, or within portable devices that have limited 
computational power. 

Summary of the Present Invention 

In a first broad form the present invention provides a method of determining parameter values used in 
25 impedance analysis of a subject, the method including, in a processing system: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 

b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 
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d) comparing the theoretical impedance values to the measured impedance values; 

e) determining at least one modification direction for at least one of the parameter value 
estimates in accordance with the results of the comparison; and, 

f) modifying at least one parameter value estimate in accordance with the determined direction, 
5 the parameter values being determined at least in part from one or more modified parameter 

value estimates. 

T3T)ically the method includes, in the computer system, iteratively modifying the parameter value 
estimates to thereby determine the parameter values. 

Typically the method includes, in the processing system, comparing the theoretical impedance values 
10 to the measured impedance values : 

a) determining a theoretical impedance locus using the parameter value estimates; 

b) determining a measured impedance locus associated with each measured impedance; and, 

c) determining an intersection between each measured impedance locus and the theoretical 
impedance locus. 

1 5 Typically the method includes, in the processing system: 

a) determining a point associated with each intersection; 

b) determining a linear regression based at least partially on the points; and, 

c) determining the modification direction based at least partially on at least one of: 
i) the relative position of the regression line and at least one of the points; and, 

20 ii) the gradient of the regression line. 

Typically the method includes, in the computer system: 

a) determining a pair of chords associated with each point of intersection, the chords 
representing the distance to locations at which the theoretical impedance locus represents 
predetermined impedance values; and, 
25 b) determining a point associated with each pair of chords; and, 

c) determining the linear regression in accordance with the value of the points. 

Typically the method includes, in the processing system: 

a) determining four impedance measurements, each impedance measurement being made at a 
respective frequency; 

30 b) comparing each measured impedance value to a corresponding theoretical impedance value; 

c) determining the modification direction based at least in part on the results of the comparison. 
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Typically the method includes, in the processing system: 

a) detennining a point based on the comparison of each measured impedance value with a 
corresponding theoretical impedance value; 

b) determining a linear regression based at least partially on the points; and, 

5 c) determining the modification direction based at least partially on at least one of: 

i) the relative position of the regression line and at least one of the points; and, 

ii) the gradient of the regression line. 

Typically the method includes, in the processing system: 

a) determining the relative position of the regression line and a predetermined point; and, 
10 b) determining the modification direction for modifying a predetermined one of the parameter 

value estimates using the relative position. 

Typically the method based on the relative position of the regression line and at least one of the 
points. 

a) comparing the gradient of the regression line and a line based on at least one of the parameter 
1 5 value estimates; and, 

b) determining the modification direction for the at least one of the parameter value estimate 
using the results of the comparison. 

Typically the method includes, in the processing system, at least one of: 

a) modifying an estimated value of Ro using a point based on the impedance measured at a first 
20 frequency; 

b) modifying an estimated value of Roo using a point based on the impedance measured at a 
second frequency; and, 

c) modifying an estimated value of a using a gradient of the regression line, where: 

a has a value between 0 and 1 ; 
25 Ro is the theoretical impedance at zero applied angular frequency; and, 

Rx, is the theoretical impedance at infinite applied angular frequency. 

Typically the method includes, in the processing system, estimating values for parameters Ro, Rx,, and 
a, where: 

a has a value between 0 and 1; 
30 Ro is the theoretical impedance at zero applied angular fi-equency; and, 

Rco is the theoretical impedance at infinite applied angular frequency. 



wo 2007/014417 



PCT/AU2006/001057 



-4- 



Typically the method includes, in the processing system, determining the theoretical impedance using 

the equation: 

_ „ Ro — i2oo 

Z = icooH where: 

Z is the measured impedance at angular frequency oo, 
5 Ro is the resistance at zero frequency, 

Rco is the resistance at infinite frequency, 

T is a time constant, and 

a has a value between 0 and 1 . 

Typically the method includes, in the processing system, using the parameter values to determine at 
10 least one of: 

a) a ratio of extra- to intra- cellular fluid; and, 

b) an indication of the at least one of the presence, absence or degree of tissue oedema in the 
subject. 

Typically the method includes, in the computer system: 
15 a) determining parameter values associated with first and second body segments; 

b) for each body segment, determining an index indicative of a ratio of the extra- to intra- 
cellular fluid; 

c) determining an index ratio based on the index for the first and second body segtnents; . 

d) comparing the index ratio to at least one reference; and, 

20 e) determining the presence, absence or degree of tissue oedema based On the results of the 

comparison. 

Typically the method includes, in the computer system, displaying an indication of at least one of: 

a) the parameter values; 

b) a ratio of extra- to intra- cellular fluid; and, 

25 c) an indication of the at least one of the presence, absence or degree of tissue oedema in the 

subject. 

Typically the method includes, in the processing system, receiving the impedance measurements from 
a measuring device, ' 

Typically the method includes, in the processing system: 
30 a) causing one or more electrical signals to be applied to the subject using a first set of 

electrodes, the one or more electrical signals having a plurality of frequencies; 
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b) determining an indication of electrical signals measured across a second set of electrodes 
applied to the subject in response to the applied one or more signals; 

c) determining from the indication and the one or more applied signals, an instantaneous 
impedance value at each of the plurality of frequencies; and, 

5 d) determining the index using the instantaneous impedance values. 

Typically the method includes, in the processing system: 

a) generating control signals, the control signals being used to apply one or more signals to the 
subject; 

b) receiving an indication of the one or more signals applied to the subject; 

10 c) receiving an indication of one or more signals measiu*ed across the subject; and, 

d) determining the impedance measurements using the received indications. 

In a second broad form the present invention provides apparatus for determining parameter values 
used in impedance analysis of a subject, the apparatus including a processing system for: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 
15 b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 

d) comparing the theoretical impedance values to the measured impedance values; 

e) determining at least one modification direction for at least one of the parameter value 
estimates in accordance with the results of the comparison; and, 

20 f) modifying the at least one parameter value estimate in accordance with the determined 

direction. 

Typically the apparatus includes: 

a) a current supply for generating an alternating current at each of the number of frequencies; 

b) at least two supply electrodes for applying the generated alternating current to a subject; 
25 c) at least two measurement electrodes for detecting a voltage across the subject; and, 

d) a sensor coupled to the measurement electrodes for determining the voltage, the sensor being 
coupled to the processing system to thereby allow the processing system to determine the 
measured impedance. 

Typically the current supply generates the alternating current at each frequency by at least one of: 
30 a) superposing a number of signals to thereby generate an alternating at each frequency 

simultaneously; and. 
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b) generating a number of alternating currents, each alternating current being at a respective 
frequency, in turn. 

Typically the apparatus is for performing the method of the first broad form of the invention. 

In a third broad form the present invention provides a method of diagnosing the presence, absence or 
5 degree of one or more conditions in a subject, the method including, in a processing system: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 

b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 

d) comparing the theoretical impedance values to the measured impedance values; 

10 e) determining at least one modification direction for at least one of the parameter value 

estimates in accordance with the results of the comparison; 
f) modifying at least one parameter value estimate in accordance with the determined direction, 
the parameter values being determined at least in part from one or more modified parameter 
value estimates; and, 

15 g) using the determined parameter values in diagnosing the presence, absence or degree of one 

or more conditions in the subject. 

It .will be appreciated that the broad forms of the invention may be used individually or in 
combination, and may be used for diagnosis of the presence, absence or degree of a range of 
conditions and illnesses, including, but not limited to oedema, pulmonary oedema, lymphodema, body 
20 composition, cardiac function, and the like. 

Brief Description of the Drawings 

An example of the present invention will now be described with reference to the accompanying 
drawings, in which: - 

Figure 1 is a schematic of an example of a theoretical equivalent circuit for biological tissue; 
25 Figure 2 is an example of a locus of impedance known as a "Wessel" plot; 

Figure 3 is a flow chart of an overview of a process for determining parameters relating to impedance 
measurements; 

Figure 4 is a schematic diagram of a first example of apparatus for monitoring bioimpedance; 
Figure 5 is a flowchart of an example of a process for performing impedance determination; 
30 Figure 6 is a schematic diagram of a second example of apparatus for monitoring bioimpedance; 

Figures 7A to 7C are a flow chart of an example of the process for determining parameters relating to 
impedance measurements; 
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Figure 8A is an example of the intersection of loci of measured impedance against and a theoretical 
impedance locus; 

Figure 8B is an example of the chords defined by the intersection of locus of measured impedance 
against a loci of theoretical locus; and, 
5 Figure 8C is an example of a plot to determine a regression line. 

Detailed Description of the Preferred Embodiments 

Figure 1 is an example of an equivalent circuit that effectively models the electrical behaviour of 
biological tissue. The equivalent circuit has two branches that represent current flow through 
extracellular fluid and intracellular fluid. The extracellular component of biological impedance is 
10 represented by and the intracellular component is represented by Rj. Capacitance of the cell 
membrane in the intracellular path is represented by C. 

The relative magnitudes of the extracellular and intracellular components of impedance of an 
alternating current (AC) are frequency dependent. At zero frequency the capacitor acts as a perfect 
insulator and all current flows through the extracellular fluid, hence the resistance at zero frequency, 
15 Ro, equals R^. At infinite frequency the capacitor acts as a perfect conductor and the current passes 
through the parallel resistive combination. The resistance at infinite frequency is given by Roo = 
RifU/(Ri+Re). 

Accordingly, the impedance of the equivalent circuit of Figure 1 at an angular frequency co, where 
a)=27c*frequency, is given by: 

20 Z = R^-^ ^""^^ (1) 

where: 

Rx>= impedance at infinite applied frequency = RiRe/(Ri+Re), 
Ro= impedance at zero applied frequency = R© and, 
T is the time constant of iJie capacitive circuit. 

25 However, the above represents an idealised situation which does not take into account the fact that the 
cell membrane is an imperfect capacitor. The effect of this can be modelled through the use of an 
alternative equivalent circuit that incorporates a constant phase element, often represented as a 
frequency dependent resistor and capacitor in parallel. Such a circuit has been previously described 
in the literature, for example in S. Grimnes and O. G. Martinsen "Bioimpedance and Bioelectrical 
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Basics", Academic Press London, 2000 and will not therefore be discussed in further detail. 
However, taking this into account leads to a modified model in which: 

Z = R„+-I1ZJ^ (2) 

where a has a value between 0 and 1 and can be thought of as an indicator of the deviation of a real 
5 system from the ideal model. 

It will be appreciated by persons skilled in the art that the following process can be performed in 
accordance with either of the above described models, but that the second model represented by 
equation (2) is preferably used as this more accurately reflects the impedance response of the subject. 

As explained above, the prior art approach to determining the desired values of Ro and R«> has been to 
10 make impedance measurements at multiple frequencies and to construct a section of a "Wessel" (also 
known as a "Cole-Cole") plot. The plot can be extrapolated to determine Ro, Roo, as shovm in Figure 2. 

An outline of one example for determining values of Ro and Roo for a given body segment will now be 
described with reference to Figure 3. 

In this example, at step 100, the impedance Z of a subject is measiu*ed at a number of different 
15 frequencies. At step 110, the measured impedances are used to provide initial estimates of parameter 
values, such as values for the parameters Rq, Roo, and a. This may be achieved using any one of a 
number of techniques as will be described in more detail below. 

In any event, at step 120, a degree of correlation between the estimated values and the impedance 
measurements is determined, for example, by using the estimated values to calculate theoretical 
20 impedance values Zt for each of the measurement frequencies ©i, using equation (2), These values 
can then be compared to the measured values of impedance Zi. 

At step 130, the results of the comparison, and in particular, differences in the measured and 
theoretical impedances Zj, Zt are used to determine a direction of modification for a selected one of 
the parameters. The parameter is modified in accordance with the determined direction, with the 
25 theoretical impedances Zt being recalculated at step 140. 



The comparison is repeated at step 150, with an assessment of whether the estimated parameter values 
are acceptable at step 160. This will ts^pically be achieved by determining if the difference between 
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the calculated and measured impedance values fall within a predetermined ranges, which is typically 
established on the basis of the level of accuracy required in the final measurement. 

In the event that the estimated values are not acceptable, the process returns to step 130, to determine 
a modification direction for either the same or a different one of the parameters. 

5 This process is repeated until the estimated parameter values are deemed to be acceptable at step 160, 
at which point the estimated values of Ro and R«, can be used in determining parameters indicative of 
various body functions, such as oedema, fat-free mass, cardiac function, or the like. 

An example of apparatus suitable for performing an analysis of a subject's impedance and 
determining impedance parameter values will now be described with reference to Figure 4. 

10 As shown the apparatus includes a measuring device 1 including a processing system 2 coupled to a 
signal generator 1 1 and a sensor 12. hi use the signal generator 1 1 and the sensor 12 are coupled to 
respective electrodes 13, 14, 15, 16, provided on a subject S, via leads L, as shown. An optional 
external interface 23 can be used to couple the measuring device 1 to one or more peripheral devices 
4, such as an external database, remote processing system or computer system, barcode scanner, or 

15 the like. 

In use, the processing system 2 is adapted to generate control signals, which cause the signal 
generator 11 to generate one or more alternating signals, such as voltage or current signals, which can 
be applied to a subject S, via the electrodes 13, 14, The sensor 12 then determines the voltage across 
or current through the subject S using the electrodes 15, 16 and transfers appropriate signals to the 
20 processing system 2. 

Accordingly, it will be appreciated that the processing system 2 may be any form of processing 
system which is suitable for generating appropriate control signals and interpreting an indication of 
measured signals to thereby determine the subject's bioelectrical impedance, and optionally determine 
other information such as cardiac parameters, or the presence absence or degree of pulmonary 
25 oedema. 

The processing system 2 may therefore be a suitably programmed computer system, such as a laptop, 
desktop, PDA, smart phone or the like. Alternatively the processing system 2 may be formed from 
specialised hardware. Similarly, the I/O device may be of any suitable form such as a touch screen, a 
keypad and display, or the like. 
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It will be appreciated that the processing system 2, the signal generator 1 1 and the sensor 12 may be 
integrated into a common housing and therefore form an integrated device. Alternatively, the 
processing system 2 may be coimected to the signal generator 11 and the sensor 12 via wired or 
wireless connections. This allows the processing system 2 to be provided remotely to the signal 
5 generator 1 1 and the sensor 12. Thus, the signal generator 1 1 and the sensor 12 may be provided in a 
unit near, or worn by the subject S, whilst the processing system 12 is situated remotely to the subject 
S. 

In one example, the outer pair of electrodes 13, 14 are placed on the thoracic and neck region of the 
subject S. However, this depends on the nature of the analysis being performed. Thus, for example, 
10 whilst this electrode arrangement is suitable for cardiac function analysis, in lymphoedema, the 
electrodes would typically be positioned on the limbs, as required. 

Once the electrodes are positioned, an alternating signal is applied to the subject S. This may be 
performed either by applying an alternating signal at a plurality of jfrequencies simultaneously, or by 
applying a number of alternating signals at different frequencies sequentially. The frequency range of 
1 5 the applied signals may also depend on the analysis being performed. 

In one example, the applied signal is a frequency rich current from a current source clamped, or 
otherwise limited, so it does not exceed the maximum allowable subject auxiliary current. However, 
alternatively, voltage signals may be applied, with a current induced in the subject being measured. 
The signal can either be constant current, impulse function or a constant voltage signal where the 
20 current is measured so it does not exceed the maximum allowable subject auxiliary current. 

A potential difference and/or current are measured between an, inner pair of electrodes 15, 16. The 
acquired signal and the measured signal will be a superposition of potentials generated by the human 
body, such as the ECG, and potentials generated by the applied current. 

Optionally the distance between the inner pair of electrodes may be measured and recorded. 
25 Similarly, other parameters relating to the subject may be recorded, such as the height, weight, age, 
sex, health status, any interventions and the date and time on which they occurred. Other information, 
such as current medication, may also be recorded. 

To assist accurate measurement of the impedance, buffer circuits may be placed in connectors that are 
used to connect the voltage sensing electrodes 15, 16 to the leads L. This ensures accurate sensing of 
30 the voltage response of the subject S, and in particular helps eliminate contributions to the measured 
voltage due to the response of the leads L, and reduces signal loss. 
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This in turn greatly reduces artefacts caused by movement of the leads L, which is particularly 
important during dialysis as sessions usually last for several hours and the subject will move around 
and change positions during this time. 

A further option is for the voltage to be measured differentially, meaning that the sensor used to 
5 measure the potential at each electrode 15 only needs to measure half of the potential as compared to 
a single ended system. 

The current measurement system may also have buffers placed in the coimectors between the 
electrodes 13,14 and the leads L. In one example, current can also be driven or sourced through the 
subject S symmetrically, which again greatly reduced the parasitic capacitances by halving the 
10 common-mode current. Another particular advantage of using a symmetrical system is that the 
micro-electronics built into the connectors for each electrode 13, 14 also removes parasitic 
capacitances that arise when the subject S, and hence the leads L move. 

The acquired signal is demodulated to obtain the impedance of the system at the applied frequencies. 
One suitable method for demodulation of superposed frequencies is to use a Fast Fourier Transform 
15 (FFT) algorithm to transform the time domain data to the frequency domain. This is typically used 
when the applied current signal is a superposition of applied frequencies. Another technique not 
requiring windowing of the measured signal is a sliding window FFT. 

In the event that the applied current signals are. formed from a sweep of different frequencies, then it 
is more typical to use a processing technique such as multiplying the measured signal with a reference 
20 sine wave and cosine wave derived from the signal generator, or with measured sine and cosine 
waves, and integrating over a whole number of cycles. This process rejects any harmonic responses 
and significantly reduces random noise. 

Other suitable digital and analog demodulation techniques will be known to persons skilled in the 
field. 

25 Impedance or admittance measurements are determined from the signals at each frequency by 
comparing the recorded voltage and current signal. The demodulation algorithm will produce an 
amplitude and phase signal at each frequency. 

An example of the operation of the apparatus for performing bioimpedance analysis will now be 
described with reference to Figure 5. 



wo 2007/014417 



PCT/AU2006/001057 



-12- 

At step 500, the processing system 2 operates to generate control signals which are provided to the 
signal generator 11 at step 510, thereby causing the signal generator to apply an alternating current 
signal to the subject S, at step 520. Typically the signal is applied at each of a number of frequencies 
fi to allow multiple frequency analysis to be performed. 

5 At step 530 the sensor 12 senses voltage signals across the subject S. At step 540 the measuring 
device, operates to digitise and sample the voltage and current signals across the subject S, allowing 
these to be used to determine instantaneous bioimpedance values for the subject S at step5. 

A specific example of the apparatus will now be described in more detail with respect to Figure 6. 

In this example, the processing system 2 includes a first processing system 10 having a processor 20, 
10 a memory 21, an input/output (I/O) device 22, and an external interface 23, coupled together via a bus 
24. The processing system 2 also includes a second processing system 17, in the form of a processing 
module. A controller 19, such as a micrologic controller, may also be provided to control activation 
of the first and second processing systems 10, 17. 

In use, the first processing system 10 controls the operation of the second processing system 17 to 
15 allow different impedance measurement procedures to be implemented, whilst the second processing 
system 17 performs specific processing tasks, to thereby reduce processing requirements on the first 
processing system 10. 

Thus, the generation of the control signals, as well as the processing to determine instantaneous 
impedance values is performed by the second processing system 17, which may therefore be formed 
20 from custom hardware, or the. like. In one particular example, the second processing system 17 is 
formed from a Field Programmable Gate Array (FPGA), although any suitable processing module, 
such as a magnetologic module, may be used. 

The operation of the first and second processing systems 10, 17, and tlie controller 19 is typically 
controlled using one or more sets of appropriate instructions. These could be in any suitable form, 
25 and may therefore include, software, firmware, embedded systems, or the like. 

The controller 19 typically operates to detect activation of the measuring device through the use of an 
on/off switch (not shown). Once the controller detects device activation, the controller 19 executes 
predefined instructions, which in turn causes activation of the first and second processing systems 10, 
17, including controlling the supply of power to the processing systems as required. 
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The first processing system 10 can then operate to control the instructions, such as the firmware, 
implemented by the second processing system 17, which in turn alters the operation of the second 
processing system 17. Additionally, the first processing system 10 can operate to analyse impedance 
determined by the second processing system 17, to allow biological parameters to be determined, 
5 Accordingly, the first processing system 10 may be formed from custom hardware or the like, 
executing appropriate applications software to allow the processes described in more detail below to 
be implemented. 

It will be appreciated that this division of processing between the first processing system 10, and the 
second processing system 17, is not essential, but there are a number of benefits that will become 
1 0 apparent from the remaining description. 

In this example, the second processing system 17 includes a PCI bridge 31 coupled to programmable 
module 36 and a bus 35, as shown. The bus 35 is in turn coupled to processing modules 32, 33, 34, 
which interface with ADCs (Analogue to Digital Converters) 37, 38, and a DAC (Digital to Analogue 
Converter) 39, respectively. 

15 The programmable module 36 is formed from programmable hardware, the operation of which is 
controlled using the instructions, which are typically downloaded from the first processing system 10. 
The firmware that specifies the configuration of hardware 36 may reside in flash memoiy (not 
shown), in the memory 21, or may be downloaded from an external source via the external interface 
23. 

• 20 Alternatively, the instructions may be stored within inbuilt memory on the second processing system 
17. In this example, the fu*st processing system 10 typically selects firmware for implementation, 
before causing this to be implemented by the second processing system 17. This may be achieved to 
allow selective activation of ftmctions encoded within the firmware, and can be performed for 
example using configuration data, such as a configuration file, or instructions representing 
25 applications software or fumware, or the like, as will be described in more detail below. 

In either case, this allows the first processing system 1 0 to be used to control operation of the second 
processing system 17 to allow predetermined current sequences to be applied to the subject S. Thus, 
for example, different firmware would be utilised if the current signal is to be used to analyse the 
impedance at a number of frequencies simultaneously, for example, by using a current signal formed 
30 from a number of superposed frequencies, as compared to the use of current signals appHed at 
different frequencies sequentially. 



wo 2007/014417 



PCT/AU2006/001057 



- 14- 

This allows a range of different current sequences to be applied to the subject by making an 
appropriate measurement type selection. Once this has been performed, the FPGA operates to 
generate a sequence of appropriate control signals I+, which are applied to the subject S. The 
voltage induced across the subject being sensed using the sensor 12, allowing the impedance values to 
5 be determined and analysed by the second processing system 17. 

Using the second processing system 17 allows the majority of processing to be performed using 
custom configured hardware. This has a number of benefits. 

Firstly, the use of a second processing system 17 allows the custom hardware configuration to be 
adapted through the use of appropriate firmware. This in turn allows a single measuring device to be 
1 0 used to perform a range of different types of analysis. 

Secondly, this vastly reduces the processing requirements on the first processing system 10. This in 
turn allows the first processing system 10 to be implemented using relatively straightforward 
hardware, whilst still allowing the measuring device to perform sufficient analysis to provide 
interpretation of the impedance. This can include for example generating a **Wesser* plot, using the 
15 impedance values to determine parameters relating to cardiac fimction, as well as determining the 
presence or absence of pulmonary oedema. 

Thirdly, this allows the measuring device 1 to be updated. Thus for example, if an improved analysis 
algorithm is created, or an improved current sequence determined for a specific impedance 
measurement type, the measuring device can be updated by downloading new firrnware via flash 
20 memory (not shovvn) or the external interface 23 . 

It will be appreciated that in the above examples, the processing is performed partially by the second 
processing system 17, and partially by the first processing system 10. However, it is also possible for 
processing to be performed by a single element, such as an FPGA, or a more generalised processing 
system. 

25 As the FPGA is a customisable processing system, it tends to be more efficient in operation than a 
more generic processing system. As a result, if an FPGA alone is used, it is generally possible to use a 
reduced overall amount of processing, allowing for a reduction in power consumption and size. 
However, the degree of flexibility, and in particular, the range of processing and analysis of the 
impedance which can be performed is limited. 
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Conversely, if only a generic processing system is used, the flexibility is enhanced at the expensive of 
a decrease in efficiency, and a consequent increase in size and power consumption. 

An example of the process for analysing impedance measurements utilising the apparatus of Figures 4 
or 6, will now be described with reference to Figures 7A to 7C, and 8Ato 8C. 

5 In particular, at step 200 the process involves measuring the impedance Zi at each of at least four 
angular frequencies a)i. At step 210, values of Ro, R«> and a are estimated and then used to calculate a 
theoretical impedance locus Zt, based on the estimated values and using equation (2) above, at step 
220. 

The initial estimation of the parameter values could be achieved in any one of a number of ways, and 
10 may include for example using a look-up table which stores estimated parameter values for given 
ranges of impedance measurements. Thus, for example, the look-up table could indicate for a given 
measurement frequency, a number of different ranges of measured impedance values, with a 
respective set of estimated values being assigned for each range. 

Alternatively, an operator may provide the initial parameter value estimates, or alternatively, pre-set 
1 5 initial estimate values may be used, depending on the preferred implementation. 

It should be noted that the theoretical impedance locus can be calculated without requiring a value for 
the parameter t on the basis that the impedance locus is a semi-circle arranged as shown in Figure 8B. 
• Thus, given values for the parameters Ro, Rco and a, then a value for x can be derived using known 
properties of equation (2) above. 

20 Given that the above described process determines values for four unknown variables, namely Ro, Roo, 
a, and x, it is necessary to use at least four measured impedance values in order to determine a unique 
solution. In this case, as the value x can be derived from the estimated values of Ro, Rco, a, the value 
of X can be used to assess the accuracy of the estimates for Ro, Roo, a. 

In this regard, the value of x is related to the characteristic frequency (the frequency of maximum 
25 impedance response) by = l/27tx. Accordingly, once values of Ro, R « and a, have been 
determined, a separate value could be determined for the characteristic frequency and hence x from 
each measured impedance value. If the parameter values for Ro, R « and a, are correct, the values of 
the characteristic frequency calculated for each measured impedance signal would be the same, 
thereby allowing the estimated values to be checked. In practice, this is achieved by creating a plot as 
30 will now be described. 
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In particular, the intersection of each impedance locus Zj with the theoretical impedance locus Zt is 
determined at step 230. An example of this is shown in Figure 8 A in which a locus for each of the 
four measured impedance values Zi, Z2, Za, Z4 is shown to intersect with the theoretical impedance 
locus Zt, At step 240, this is used to determine chord lengths Ui, Vi for each intersection. An example 
5 of how the chord lengths Uj, Vi are determined is shown in Figure 8B. 

It will be appreciated by persons skilled in the art that whilst steps 220 to 240 could be achieved by 
generating a physical plot of the measured and theoretical impedance loci Zi, Zt, this is not essential. 
Thus, for example, suitable calculations can alternatively be performed to determine the theoretical 
coordinates of the intersections. 

10 At step 250, a chord length ratio UiA^i is deteniiined for each intersection, with this being used to plot 
corresponding points Pi, P2, Pa, P4 by plotting loge(UiA^i) loge((»i), for each of the intersections, at step 
260. At step 270 the processing system 2 performs a linear regression utilising the points Pj, P2, P3, 
P4 to determine a regression line L shown in Figure SC. 

Whist a regression line L is used in this example, the process could operate using any line extending 
15 between Pj and P4 as the technique examines the gradient of the line and then modifies this until it fits 
with observed values. However, the use of a regression line ensures that the initial fit is as accurate as 
possible, thereby reducing computational requirements. 

The processing system 2 can then use the regression line to improve the estimated values of Ro, R<» 
and a. This can be achieved taking into account that: 
20 • The position of point P2 above or below the regression line is sensitive mainly to changes in 

Ro, and can therefore be used as feedback for Rq estimation. 

• The position of point P3 above or below the regression line is sensitive mainly to changes in 
Roo, and can therefore be used as feedback for R«, estimation. 

• The slope of the regression line L can be compared with a target slope T defined as (a-1), and 
25 can therefore be use to provide feedback for a estimation. 

In this regard, the value of points P2 and P3 are used for estimation of Ro and R«, respectively as the 
points are based on the impedance measurements near to Ro and Roo, and therefore have a greater 
influence on these values. 



At step 280 the processing system 2 determines the relative position of the point P2 relative to the 
30 regression line L and determines if the point P2 lies on the regression line L at step 290. 
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In this regard, the assessment will typically not be absolute, and will take into account a required 
precision for the resulting values of Ro, Roo. Thus, the assessment will typically determine if the point 
P2 is within a predetermined distance from the regression line i. 

If not, the processing system 2 operates to revise the value of Ro and recalculates the regression line at 
5 step 300. The revision of Ro is performed so that whether the value of Ro is increased or decreased 
depends on the relative position of the point P2 to the regression line £, Thus, for example, if P2 is 
above the line then Ro is decreased in value. This ensures that the modification of the estimated value 
of Ro will always result in an improved estimate, thereby reducing the mathematical requirements for 
calculating suitable estimates. 

10 The recalculation of the regression line L is performed using the general method set out in steps 210 
to 270 and this will not therefore be described in detail. 

If the point P2 lies on the regression line Z, or once the regression line has been revised at step 300, 
the processing system 2 determines the position of the point P3 relative to the regression line L at step 
310. Again if it is determined that the point P3 does not lie on, or is not within a predetermined 
15 distance of, the regression line L at step 320, this is used to revise the value or R«, and recalculate the 
regression line at step 330. In this case, if P3 is above the line then R co is decreased in value. Again 
this is achieved using the steps by repeating the steps 210 to 270. 

Once this has been completed, the processing system 2 assesses whether all four points Pi, P25 P3, P4 
lie on, or within a predetermined distance of the regression line Z, at step 340. If this is not the case, 
20 the processing system 2 will return to step 280 to determine the position of point P2 and repeat the 
steps 280 to 340. 

Once all four points lie on, or within a predetermined distance of the regression line Z, the process 
moves on to step 350 allowing the processing system 2 to compare the gradient of the regression line 
L to the gradient of the target slope T, which is given by a - i. 

25 If it is assessed that the gradients are not equal, or within a predetermined range, at step 360, the 
processing system 2 revises the value of a and returns to step 280 to recalculate the regression line Z, 
and again assess the relative position of the points P2, P3 and the regression line Z. 

Once values for the estimates of Ro, R «> and a are finalised, the points Pi, P2, P3, P4 should lie on a 
straight line with a slope of a-1, with log QJIW) being zero at the characteristic frequency/: = 1/2 nr. 
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Thus, by using the estimated values to determine a common regression line Z, this effectively tests the 
accuracy of the estimates for Ro, R«, and a. 

At step 380 the processing system 2 determines at least the values of Ro and Ro, and can optionally 
also determine the value x. These values can then be used in assessing the subject. 

5 This can be achieved in a number of manners depending on the particular aspect of the patient being 
monitored. 

For example, in assessing the presence, absence or degree of tissue oedema, it is typical to use the 
parameters to determine an index I based on a ratio of extra- to intra- cellular fluid, given by the 
equation: 

10 = (3)- 

This is typically achieved by comparing the index / to a predetermined reference or range, v^hich may 
be derived for example from: 

• study of a reference population; 

• comparison to an index previously determined for the subject, known as a longitudinal 
15 analysis; and, 

• an index obtained from one or more other body segments for the individual. 

Thus, the parameter values derived using the above techniques may be used in studying body 
composition, oedema detection, monitoring cardiac function, or the like. 

hi any event, the above-described process operates by estimating parameter values that can be used to 
20 model the impedance. The theoretical impedance values predicted by the model are compared to the 
measured values, allowing the estimated parameter values to be revised, hi this regard, the 
comparison is performed in such a manner as to determine the direction in which the estimated 
parameter values should be revised, thereby vastly reducing the processing required to achieve 
acceptable parameter values. 

25 Thus, this should be contrasted with other iterative techniques, such as method of least squares 
approaches, which only assess the magnitude of differences between calculated and measured values. 
In such approaches, assessment of magnitude only means that modification of the estimated values 
can result in worse estimates, thereby increasing the processing required to achieve a satisfactory 
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estimate. Furthermore, because, the process effectively adjusts each of four parameter values 
corresponding to Ro, Ra, x and a, and uses this to "explore" an error surface, there is a risk that the 
process will locate a local minimum value that does not necessarily corresponding to the best 
parameter combination. 

5 In contrast, the above process derives feedback for controlling the adjustment direction of three 
estimated parameter values Ro, Roo, and a, with the value of the fourth parameter value t being used to 
ensure the estimated values are correct. Thus, as direction of modification of the estimated values is 
determined, this vastly reduces processing requirements, making the process suitable for determining 
the values Ro, Roo, in real-time. 

10 In the above described example, the processing of the measured impedance values is described as 
being performed by the processing system 2. In this regard, it will be appreciated that the process may 
be performed entirely by either one of the first or second processing systems 10, 17, or a combination 
or the two. 

Additionally, or alternatively, any one or more of processing tasks may be performed by a remote 
15 processing system, such as a computer system or the like, coupled to the measuring device via the 
external interface 23. In this arrangement, the remote processing system therefore effectively forms 
part or all of the processing system 2. Accordingly, this allows the measuring device 1 to perform the 
impedance measurement procedure and provide the resulting impedance measurements to the remote 
processing system, to allow the impedance parameter values to be calculated using the above 
20 described process. Alternatively, calculation may be performed entirely within the measuring device. 

In any event, it will be appreciated that in one example, a remote processing system may perform 
tasks otherwise performed by the first processing system 10, allowing the measuring device 1 to be 
provided including only the second processing system 17 coupled directly to the external interface 23. 
This allows the measuring device 1 to be controlled by the remote processmg system. This would 
25 typically be achieved via the use of suitable applications software installed on the remote processing 
system. 

In another example, the above described process may be used with an alternative measuring devices, 
including implantable devices such as cardiac pacemakers or defibrillators. In particular, in the case 
of implantable devices, it is typical to minimise both the power consumption and physical size of the 
30 device, to make the device practical for use. As a result, implantable devices typically have only 
limited processing ability, and are not generally able to measure the phase of signals at more than one 
frequency. 
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Accordingly, the above described process could be used to estimate various cardiac parameters such 
as oedema or heart function by measuring the magnitude of impedance only at several frequencies. 
This approach would enable the estimation of multifrequency results for low power small implantable 
devices, 

5 The discussion has referred to both oedema and lymphoedema, as it is clear to a skilled addressee that 
the above method and apparatus may be utilised on any form of impedance measurements including, 
but not limited to body composition, cardiac function, tissue oedema, pulmonary oedema, or the like. 
The method may also be used to derive indices that can be compared to references, such as previously 
determined indices, indices derived from different anatomical region or different individuals, thereby 
10 allowing different forms of analysis to be performed, depending on the circumstances in, which the 
process is used. 

Throughout the specification, the aim has been to describe the preferred embodiments of the invention 
without limiting the invention to any one embodiment or specific collection of features. Various 
changes arid modifications may be made to the embodiments described and illustrated without 
1 5 departing firom the present invention. 

Thus, for example, it will be appreciated that features from different examples above may be used 
interchangeably where appropriate. Furthermore, whilst the above examples have focussed on a 
subject such as a huinan, it will be appreciated that the measuring device and techniques described 
above can be used with any animal, including but not limited to, primates, livestock, performance 
20 animals, such as race horses, or the like. 

The above described processes can be used for diagnosing the presence, absence or degree of a range 
of conditions and illnesses, including, but not limited to oedema, lymphodema, body composition, 
cardiac function, or the like. 
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THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS : 

1) A method of determining parameter values used in impedance analysis of a subject, the method 
including, in a processing system: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 
5 b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 

d) comparing the theoretical impedance values to the measured impedance values; 

e) determining at least one modification direction for at least one of the parameter value 
estimates in accordance w^ith the results of the comparison; and, 

10 f) modifying at least one parameter value estimate in accordance w^ith the determined direction, 

the parameter values being determined at least in part from one or more modified parameter 
value estimates. 

2) A method according to claim 1, wherein the method includes, in the computer system, iteratively 
modifying the parameter value estimates to thereby determine the parameter values. 

15 3) A method according to claim 1 or claim 2, wherein the method includes, in the processing system, 
comparing the theoretical impedance values to the measured impedance values: 

a) determining a theoretical impedance locus using the parameter value estimates; 

b) determining a measured impedance locus associated with each measured impedance; and, 

c) determining an intersection between each measured impedance locus and the theoretical 
20 impedance locus. 

4) A method according to claim 3, wherein the method includes, in the processing system: 

a) determining a point associated with each intersection; 

b) determining a linear regression based at least partially on the points; and, 

c) determining the modification direction based at least partially on at least one of: 
25 i) the relative position of the regression line and at least one of the points; and, 

ii) the gradient of the regression line. 

5) A method according to claim 4, wherein the method includes, in the computer system: 

a) determining a pair of chords associated with each point of intersection, the chords 
representing the distance to locations at which the theoretical impedance locus represents 

30 predetermined impedance values; and, 

b) determining a point associated with each pair of chords; and, 

c) determining the linear regression in accordance with the value of the points. - 

6) A method according to any one of the claims 1 to 5 1, wherein the method includes, in the 
processing system: 
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a) determining four impedance measurements, each impedance measurement being made at a 
respective frequency; 

b) comparing each measured impedance value to a corresponding theoretical impedance value; 

c) determining the modification direction based at least in part on the results of the comparison. 
5 7) A method according to claim 6, wherein the method includes, in the processing system: 

a) determining a point based on the comparison of each measured impedance value with a 
corresponding theoretical impedance value; 

b) determining a linear regression based at least partially on the points; and, 

c) determining the modification direction based at least partially on at least one of: 
10 i) the relative position of the regression line and at least one of the points; and, 

ii) the gradient of the regression line. 

8) A method according to claim 7, wherein the method includes, in the processing system: 

a) determining the relative position of the regression line and a predetermined point; and, 

b) determining the modification direction for modifying a predetermined one of the parameter 
1 5 value estimates using the relative position. 

9) A method according to claim 7 or claim 8, wherein the method based on the relative position of 
the regression line and at least one of the points. 

a) comparing the gradient of the regression line and a line based on at least one of the parameter 
value estimates; and, 

20 b) determining the modification direction for the at least one of the parameter value estimate 

using the results of the comparison. 

10) A method according to claim 7, wherein the method includes, in the processing system, at least 
one of: 

a) modifying an estimated value of Ro using a point based on the impedance measured at a first 
25 frequency; 

b) modifying an estimated value of Roo using a point based on the impedance measured at a 
second frequency; and, 

c) modifying an estimated value of a using a gradient of the regression line, where: 

a has a value between 0 and 1; 
30 Ro is the theoretical impedance at zero applied angular frequency; and, 

Roo is the theoretical impedance at infinite applied angular frequency. 

11) A method according to any one of the claims 1 to 10, wherein the method includes, in the 
processing system, estimating values for parameters Ro, Roo, and a, where: 

a has a value between 0 and 1 ; 
35 Ro is the theoretical impedance at zero applied angular frequency; and. 
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Rco is the theoretical impedance at infinite applied angular frequency. 

12) A method according to any one of the claims 1 to 11, wherein the method includes, in the 

processing system, determining the theoretical impedance using the equation: 

^ „ Ro-Roo 

Z = Roo+ where: 

5 Z is the measured impedance at angular frequency co, 

Ro is the resistance at zero frequency, 
Roo is the resistance at infinite frequency, 
T is a time constant, and 
a has a value between 0 and 1 . 
10 13) A method according to any one of the claims 1 to 12, wherein the method includes, in the 
processing system, using the parameter values to determine at least one of: 

a) a ratio of extra- to intra- cellular fluid; and, 

b) an indication of the at least one of the presence, absence or degree of tissue oedema in the 
subject. 

15 14) A method according to any one of the claims 1 to 13, wherein the includes, in the computer 
system: 

a) determining parameter values associated with first and second body segments; 

b) for each body segment, determining an index indicative of a ratio of the extra- to intra- 
cellular fluid; 

20 c) determining an index ratio based on the index for the first and second body segments; 

d) comparing the index ratio to at least one reference; and, 

e) determining the presence, absence or degree of tissue oedema based on the results of the 
comparison. 

15) A method according to claim 13 or claim 14, wherein the method includes, in the computer 
25 system, displaying an indication of at least one of: 

a) the parameter values; 

b) a ratio of extra- to intra- cellular fluid; and, 

c) an indication of the at least one of the presence, absence or degree of tissue oedema in the 
subject. 

30 16) A method according to any one of the claims 1 to 15, wherein the method includes, in the 
processing system, receiving the impedance measurements from a measuring device. 
17) A method according to any one of the claims 1 to 16, wherein the method includes, in the 
processing system: 
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a) causing one or more electrical signals to be applied to the subject using a first set of 
electrodes, the one or more electrical signals having a plurality of frequencies; 

b) determining an indication of electrical signals measured across a second set of electrodes 
applied to the subject in response to the applied one or more signals; 

5 c) determining from the indication and the one or more applied signals, an instantaneous 

impedance value at each of the plurality of frequencies; and, 
d) determining the index using the instantaneous impedance values. 
18) A method according to any one of the claims 1 to 16, wherein the method includes, in the 
processing system: 

10 a) generating control signals, the control signals being used to apply one or more signals to the 

subject; 

b) receiving an indication of the one or more signals applied to the subject; 

c) receiving an indication of one or more signals measured across the subject; and, 

d) determining the impedance measurements using the received indications. 

15 19) Apparatus for determining parameter values used in impedance analysis of a subject, the 
apparatus including a processing system for: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 

b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 
20 d) comparing the theoretical impedance values to the measured impedance values; 

e) determining at least one modification direction for at least one of the parameter value 
estimates in accordance with the results of the comparison; and, 

f) modifying the at least one parameter value estimate in accordance with the determined 
direction. 

25 20) Apparatus according to claim 19, wherein the apparatus includes: 

a) a current supply for generating an alternating current at each of the number of frequencies; 

b) at least two supply electrodes for applying the generated alternating current to a subject; 

c) at least two measurement electrodes for detecting a voltage across the subject; and, 

d) a sensor coupled to the measurement electrodes for determining the voltage, the sensor being 
30 coupled to the processing system to thereby allow the processing system to determine the 

measured impedance. 

21) Apparatus according to claim 20, wherein the current supply generates the alternating current at 
each frequency by at least one of: 

a) superposing a number of signals to thereby generate an alternating at each frequency 
35 simultaneously; and, 
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b) generating a number of alternating currents, each alternating current being at a respective 
frequency, in turn. 

22) Apparatus according to any one of the claims 19 to 21, wherein the apparatus is for performing 
the method of any one of the claims 1 to 18. 

23) A method of diagnosing the presence, absence or degree of one or more conditions in a subject, 
the method including, in a processing system: 

a) determining a number of impedance measurements at a corresponding number of frequencies; 

b) determining estimates of the parameter values; 

c) determining theoretical impedance values based on the parameter value estimates; 

d) comparing the theoretical impedance values to the measured impedance values; 

e) determining at least one modification direction for at least one of the parameter value 
estimates in accordance with the results of the comparison; 

f) modifying at least one parameter value estimate in accordance with the determined direction, 
the parameter values being determined at least in part from one or more modified parameter 
value estimates; and, 

g) using the determined parameter values in diagnosing the presence, absence or degree of one 
or more conditions in the subject. 
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